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Abstract. We propose a unifying phase-space approach to the construction 
of mutually unbiased bases for a two-qubit system. It is based on an 
explicit classification of the geometrical structures compatible with the notion of 
unbiasedness. These consist of bundles of discrete curves intersecting only at the 
origin and satisfying certain additional properties. We also consider the feasible 
transformations between different kinds of curves and show that they correspond 
to local rotations around the Bloch-sphere principal axes. We suggest how to 
generalize the method to systems in dimensions that are powers of a prime. 



1. Introduction 

The notion of mutually unbiased bases (MUBs) emerged in the seminal work 
of Schwinger [1] and it has turned into a cornerstone of the modern quantum 
information. Indeed, MUBs play a central role in a proper understanding of 
complementarity [2-6] , as well as in approaching some relevant issues such as optimum 
state reconstruction [7,8], quantum key distribution [9, 10], quantum error correction 
codes [11,12], and the mean king problem [13-17]. 

For a d-dimensional system (also known as a qudit) it has been found that the 
maximum number of MUBs cannot be greater than d + 1 and this limit is reached 
if d = p is prime [18] or power of prime, d = p n [19]. It was shown in reference [20] 
that the construction of MUBs is closely related to the possibility of finding of d + 1 
disjoint classes, each one having d— 1 commuting operators, so that the corresponding 
eigenstates form sets of MUBs. Since then, different explicit constructions of MUBs 
in prime power dimensions have been suggested in a number of papers [21-27]. 

The phase space of of a qudit can be seen as a d x d lattice whose coordinates 
are elements of the finite Galois field GF(d) [28]. At first sight, the use of elements of 
GF(d) as coordinates could be seen as an unnecessary complication, but it proves to be 
an essential step: only by doing this we can endow the phase-space grid with the same 
geometric properties as the ordinary plane. There are several possibilities for mapping 
quantum states onto this phase space [29-31]. However, special mention must be 
made of the elegant approach developed by Wootters and coworkers in references [32] 
and [33] , which has been used to define a discrete Wigner function (see references [34] 
and [35] for picturing qubits in phase space) . Any good assignment of quantum states 
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to lines is called a 'quantum net". In fact, there is not a unique quantum net for a 
given d x d phase space. However, one can manage to construct lines and striations 
(sets of parallel lines) in this phase space: after an arbitrary choice that does not lead 
to anything fundamentally new, it turns out that the orthogonal bases associated with 
each striation are mutually unbiased. 

In this paper, we proceed just in the opposite way. We start by considering 
the geometrical structures in phase space that are compatible with the notion of 
unbiasedness. By taking the case of two identical two-dimensional systems (i.e., two 
qubits) as the thread for our approach, we classify these admissible structures into 
rays and curves (and the former also in regular and exceptional, depending on the 
degeneracy). To each bundle of curves, we associate a MUB, and we show how these 
MUBs are related by local transformations that do not change the corresponding 
entanglement properties. Finally, we sketch how to extend this theory to higher (power 
of prime) dimensions. We hope that this new method can seed light on the structure 
of MUBs and can help to resolve some of the open problems in this field. For example, 
all the MUB structures in 8- and 16-dimensional Hilbert space are known [36], but 
in the 16-dimensional case the transformations to go from one structure to any other 
are unknown and hitherto a method to find them (in any space dimension) has been 
lacking. Our approach provides a means to find such transformations in a systematic 
manner. 



2. Constructing a set of mutually unbiased bases 

When the space dimension d — p n is a power of a prime it is natural to conceive the 
system as composed of n constituents, each of dimension p [37] . We briefly summarize 
a simple construction of MUBs for this case, according to the method introduced 
in reference [27], although focusing on the particular case of two-qubits. The main 
idea consists in labeling both the states of the subsystems and the generators of the 
generalized Pauli group (acting in the four-dimensional Hilbert space) with elements 
of the finite field GF{4), instead of natural numbers. In particular, we shall denote 
as | a) with a G GF{4) an orthonormal basis in the Hilbert space of the system. 
Operationally, the elements of the basis can be labeled by powers of a primitive element 
(that is, a root of the minimal irreducible polynomial, a 2 + a + 1 = 0), so that the 
basis reads 

{|0>, \<t 2 = <t + 1), |a 3 = l)}. (2.1) 
These vectors are eigenvectors of the generalized position operators Zp 

Z = Yl X(ap)\a)(a\, (2.2) 

aeGF(4) 

where henceforth we assume a, (3 S GF(A). Here x(#) is an additive character 



X{6) = cxp 



27TI 

tr(0) 

P 



(2.3) 



and the trace operation, which maps elements of GF(4) onto the prime field GF(2) ~ 
Z 2 , is defined as tr(9) = 9 + 2 . The diagonal operators Zp are conjugated to the 
generalized momentum operators Xp 

Xp= ]T \a + p)(a\, (2.4) 
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precisely through the finite Fourier transform 

FXpF* =Z P , (2.5) 

with 

F=\ E X(a(3)\a){p\. (2.6) 

The operators {Z a: Xp} are the generators of the generalized Pauli group 

Z a X fJ = X (a/3) X p Z a . (2.7) 

In consequence, we can form five sets of commuting operators (which from now on 
will be called displacement operators) as follows, 

{ZaXf3 = ^ a } , (2.8) 

with fi <E GF(4). The displacement operators (|2.8p can be factorized into products of 
powers of single-particle operators o~ z and o~ x , whose expression in the standard basis 
of two-dimensional Hilbert space is 

<r z = |0)(0|-|l)(l|, a x = |0)(l| + |l)(0|. (2.9) 

This factorization can be carried out by mapping each element of GF(A) onto an 
ordered set of natural numbers [33], a <^ (oi, 0.2), where aj are the coefficients of the 
expansion of a in a field basis 9j 

a = a 1 6 1 +a 2 6 2 . (2.10) 

A convenient field basis is that in which the finite Fourier transform is factorized into 
a product of single-particle Fourier operators. This is the so-called self-dual basis, 
defined by the property tr(6i6j) — &y. In our case the self-dual basis is (a, a 2 ) and 
leads to the following factorizations 

Z a =a a /<j?, X fi = (2.11) 

where a = a\o~ + a 2 o 2 and f3 — b\o + b 2 <J 2 ■ Using this factorization, one can 
immediately check that, among the five MUBs that exist in this case, three are 
factorable and two are maximally entangled [38]. Although the factorization of a 
particular displacement operator depends on the choice of a basis in the field, the 
global separability properties (i.e., the number of factorable and maximally entangled 
MUBs) is basis independent. That is, any nonlocal unitary transformation that yields 
only factorable or maximally entangled bases (i.e., a transformation from the Clifford 
group) will provide an isomorphic set of MUBs with respect to the separability, except, 
perhaps, for some trivial permutations. Nevertheless, this property holds only for two 
qubits because for higher-dimensional cases more complicated structures arise [36] . 



3. Mapping the mutually unbiased bases onto phase space 

The problem of MUBs can be further clarified by an appropriate representation in 
phase space, which is defined as a collection of ordered points (a, /?) 6 GF(4) x GF(4). 
In this finite phase space the operators from the five sets (|2.8p are labeled by points 
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Table 1. Rays and their associated physical operators. 



Basis 


Ray 


Factorized operators 


1 


13 = 


ffzt 


lcr z cr z cr z 


2 


(5 = a 


a y t 


HcTy O-yCTy 


3 


f} = era 




a x (7y Uycr z 


1 


(5 = a 2 a 


CTy&x 




5 


a = 


cr x 1 


\a x a x a x 





4 


3 


2 




3 


2 


4 




2 


4 


3 




1 


1 


1 



Figure 1. Phase-space picture corresponding to the construction in table 1. 

of rays (i.e., 'straight lines passing through the origin). The vertical axis has a = 
and the horizontal axis has (3 = 0. For our case, we explicitly have 



p = 


- 




P = 


a 


"* Z a X a , Z a 2X a 2, Z a 3X a 3 


= 


act - 


-» Z a X a 2 , Z„iX a 3 , Z a zX a 


P = 


a 2 a - 


-> Z a X a 3,Z a iX a ,Z a 3X„i 


a = 


- 





where the left column indicates the ray corresponding to the operators appearing in 
the three rightmost columns. In the factorized form, the set in (|3.ip can be expressed 
as in table 1. 

In figure 1 we plot the phase-space representation of the sets of operators in 
table 1. Each set has been arbitrarily assigned to the number appearing in the left 
column of the table. The sets of operators 1 and 5 define the horizontal and the vertical 
axes, respectively, and they lead, together with the operators associated to line 2, to 
three separable bases (i.e., the three operators in each of the first three rows commute 
for each of the two subsystems, separately). In physical space, all these operators can 
be associated with rotations of each qubit around the z-, x- and y-axis, respectively. 
Eigenstates of the operators associated with the lines 3 and 4 form entangled bases (in 
fact, their simultaneous eigenstates are all maximally entangled states). The origin is 
labeled as o and is the common intersecting point of all the rays. 

It is clear that under local transformations the factorable and entangled MUBs 
preserve their separability properties. Two natural questions thus arise in this respect: 
Is the arrangement in table 1 and the corresponding geometrical association with rays 
in phase-space unique? If this is not the case, why do different arrangements always 
lead to the same separability structure of MUBs? 
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4. Curves in phase space 

To answer these questions we shall approach the problem from a different perspective, 
namely, by determining all the possible geometrical structures in phase space that 
correspond to MUBs. First of all, let us observe that any ray can be defined in the 
parametric form 

a(K)=rjK, /3(k) = ^k, (4-1) 

where tj,( e GF(A) are fixed while k 6 GF(4) is a parameter that runs through all 
the field elements. The rays (|4.1[) can be seen as the simplest nonsingular (i.e., no 
self-intersecting) Abelian substructures in phase space, in the sense that 

a( K + k') = <x{k) +ol{k') , P(k + k) = (3{k) + P(k') . (4.2) 

However, the rays are not the only Abelian structures: it is easy to see that the 
parametric curves (that obviously pass through the origin) 

a(n) — HoK + Mik 2 , (3(k) = 7] k + -qin 2 , (4.3) 

also satisfy the condition (|4.2[) . If, in addition, we impose 

tr(a/3') = tr(a'/3) , (4.4) 

where a' = a(n') and (3' = /?(«/)> then the displacement operators associated to the 
curves (|4.3|) commute with each other and the coefficients jij and rjj must satisfy the 
following restrictions (commutativity conditions) 

mm + (mm) 2 = Mom + (HoVi) 2 ■ (4-5) 

All the possible Abelian curves satisfying condition (|4.5|) can be divided into two types: 

a) regular curves 

a— curves : a = an , (3 = rjn + a 2 n 2 , 

(4.6) 

j3— curves : (3 = crn , a = tjk + ct 2 k 2 . 

b) exceptional curves 

a = fi{n + K 2 ), (3 = li 2 {(jk + a 2 K 2 ) . (4.7) 

The regular curves are nondegenerate, in the sense that a or (3 (or both) are not 
repeated in any set of four points defining a curve. In other words, a or f3 (or both) 
take all the values in the field GF(4). This allows us to write down explicit relations 
between a and (3 as follows 

a— curves : (3 = r\a 2 a + a 2 , 

(4.8) 

j3— curves : a — rja 2 (3 + (3 2 . 

By varying the parameter r\ in the first of equations (|4.8j) we can construct the a-curves 
in table 2, which show a different arrangement of operators than (|3.1[) . Figure 2 
shows the corresponding points of a-curves in phase space. Note, that we have 
completed table 2 and figure 2 with the vertical [X ai X a i , X a -s) axis. The factorization 
of operators in each table (the self-dual basis is used for the representation of operators 
in terms of Pauli matrices) is different from the standard one in table 1. The curves 
marked as 3, 4 and 5 lead now to factorable MUBs, while the ones marked as 1 and 
2 lead to maximally entangled bases. 
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Table 2. a-curves and their corresponding operators. 



Basis 


a- 


curves 


Displacement operators 


Factorized operators 


1 


a 


= a 2 


Z a 2X a , Z a 3X a 3 , Z a X a 2 




CTyUy <Jz<Jx 


2 


p 


= a + a 2 




<J X Uy 


<Jz<Jz <Jy<J x 


3 


p 


= era + a 2 


Z a 2 X a 2 , Z a 3 X a 2 , Za 


la y 


CT Z CTy <7 2 1L 


■1 




= a 2 a + a 2 


Z a 2, Z^zXa, Z a X a 


ta z 


CTyCTz 


5 


a 


= 


X a ,X a 2,X a i 


a x \ 


1L(T X a x a x 



5 


2 


2 


1 




1 


3 


3 




4 


1 


4 




3 


4 


2 



Figure 2. Phase-space picture corresponding to the construction in table 2. 



The /3-curves and the corresponding table can be obtained from table 2 by 
exchanging a and [3 (and correspondingly Z and X operators) and is given in table 3. 
The phase-space picture corresponding to table 3 is shown in figure 3 and can be easily 
obtained from figure 2 by mirroring this figure about the main diagonal. Observe that 
the curves (3 = a 2 and a = 1 then become identical, since this curve is symmetric 
about the diagonal. 



It is worth noting that all the a-curves, except (3 



are /3-degenerate: the 



same value of (3 corresponds to different values of a. Obviously, the analogous a- 
degeneration appears in the /3-curves. 

Exceptional curves (|4.7p have quite a different structure. Now, every point is 
doubly degenerate and can be obtained from equations that relate powers of a and f3: 

a 2 =(ia, (3 2 = t i 2 f3. (4.9) 



Table 3. Phase-space /3-curves and their corresponding operators. 



Basis ftefa-curves Displacement operators Factorized operators 



1 


a 


= p 2 


X a 2 Z a ,X a 3 Z a 3 , X a Z a 2 


<Jz<J x 


UyUy 


<T x Oz 


2 


a 


= p + p 2 


Jf CT 2Z cr 3,X (r 3,X CT Z CT 3 


az<?y 


<T X T X 


OyOz 


3 


a 


= aP + /3 2 


X a 2 Z <r 2 > X a3 Z a 2 > X ° 


tOy 


a x a y 


<J X 1 


4 


a 


= a 2 /3 + f3 2 


X^2 , X^Za, X a Z a 


l_cr x 


OyO X 


<Tyt 


5 


a 


= 


Z a ,Z a 2,Z a 3 


a z l 


tuz 


&Z&Z 



2 


4 


3 


1 


4 


1 


3 


2 


3 


4 


1 


2 


o 


5 


5 


5 



Figure 3. Phase-space picture corresponding to the construction in table 3. 
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Table 4. Bundle consisting of two exceptional curves, one a-curve, one /3-curve, 
and a ray. 

Basis Curves and rays Displacement operators Factorized operators 



o x a x a z a z <?yGy 

(J x \ 1<T Z cr x cr z 

IcTy (7 z (7y (T z \ 

lo-j. cr y a x (Tyl 

a z a x cr x <jy cr y u z 



1 


4 


5 


1 


4 


5 


3 


3 


2 


4 


2 


5 


O 


3 


2 


1 



Figure 4. Phase-space picture corresponding to the construction in table 4. 

It is impossible to write an explicit nontrivial equation of the form f(a,f3) = for 
them. The existence of these curves allows us to obtain interesting arrangements of 
MUB operators in tables that do not contain any axis (z, x or y). There are two of 
such structures, shown in tables 4 and 5. As can be seen from the rightmost column in 
both tables, the physical difference between the two structures is that the two qubits 
are permuted between them. The lines marked 2, 3 and 4 in both tables lead to 
factorable MUBs, while the lines marked as 1 and 5 give maximally entangled ones. 



^ a = re + re 

fj = (TR + <T 2 K 2 

a = ct 2 (k + re 2 ) 

2 /3 = CT 2 K + re 2 

3 /3 = era + a 2 

4 a = a 2 P + P 2 

5 P = era 



X a , Z a 2, Z a 2X a 

X^l , Z a X a 3 , Z a X,y 

Z a X 2 , Z iX 3 , Z sXcr 



Table 5. Bundle consisting of two exceptional curves, one a-curve, one /3-curve, 
and a ray. 



Basis Curves and rays Displacement operators Factorized operators 



1 


a = 
P = 


re + re 2 
crre + a 2 n 2 


X a -3,Z (7 3,Z a .3X a .3 


<Tx&x 


& z <y z 


CTyCTy 


2 


a = 
P = 


tr(re + re 2 ) 
re + uk 2 


X a 2 , Z a , Z CT X CT 2 


ltT x 


o z t 


a z a x 




P = 


a 2 a + a 2 


Z a 1 , Z a 3X a , ZaXcj 


ta z 


Cy<Tz 


(Tyl 


1 


a = 


aP + P 2 


z ^ x ^^ z ^X (T 3,X cy 


ta y 


<J X <Jy 


a x \ 


5 


P = 


2 

a a 


Z(jX a 3, Z a lX a , Z a 3X a 2 


<Ty(T x 


<J X <J Z 


a z a y 





5 


4 


1 


2 

k 


2 
3 
2 


4 
5 

3 


5 

3 
1 



Figure 5. Phase-space picture corresponding to the construction in table 5. 
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Table 6. Bundle consisting of two exceptional curves, one a-curve, one /3-curve, 
and a ray. 



Basis Curves and rays Displacement operators Factorized operators 



1 


a = 
P = 


<t 2 (k + re 2 ) 
a 2 re + re 2 


X^2 , Zcr, Z a X a 2 






<TzCr x 


2 


a = 
P = 


ct(k + re 2 ) 
crre + crre 2 


X a , Z a 2, Z a 2X a 


a x 1L 




<T x (Jz 


3 


f3 = 


a + a 2 




(Jy<T X 


O ~ z O ' z 


<T X (Ty 


4 


a = 


f3 + /3 2 


Z ^X iy 2,X iy 3,Z l7 3X a 






CFyCFz 


5 


f3 = 


a 


Z a X a , Z a 2X a 2 , .2^3X^3 


CTyt 


ta v 


CTyCFy 



Finally, there is a last table containing two exceptional curves and a ray 
corresponding to the spin operators in the y-direction, as it is shown in table 6. 

To sum up, there exist fifteen different Abelian structures, five rays and ten 
curves, which can be organized in six different forms with the respect to MUBs. The 
existence of only six bundles of mutually nonintersecting Abelian nonsingular curves 
(i.e., different tables) also follows from the fact that the coset of the full symplcctic 
group, which preserves the commutation relations (|2.7|) . on operations corresponding 
to nontrivial permutations of columns and rows of (any) table [generated by the 
symplectic group Sp(2, GF(4))], is precisely of order 6. 



5. The effect of local transformations 



As we have noticed, different arrangements of operators in tables (or bundling of 
phase-space curves) lead to the same separability structure. To understand this point, 
let us study the effect of local transformations. In other words, we wish to characterize 
how a given curve changes when a local transformation is applied to a set of operators 
labeled by points of this curve. 

To deal with such operations with curves, let us recall that a generic displacement 
operator is factorized in the self-dual basis as 

Z a X p - (<t>^)« 2 ^) = (aiA) ® (02,62) ■ (5.1) 

It is clear that under local transformation (rotations by tt/2 radians around the z-, 
x- or y-axes) applied to the jth particle (j — 1,2), the indices of the displacement 
operators are transformed as follows: 

z— rotation : (aj, bj) — > (aj + b^ bj) , 

x— rotation : (a.,-, bj) — * (a.,-, bj + a,j) , (5-2) 
y— rotation : (a,j, bj) — > (a,j + aj + bj, bj + a,j +bj) — (bj, ctj) . 

To give a concrete example, suppose we consider a z-axis rotation. The operator a z , 
corresponding to (aj = l,bj = 0), is transformed into (a 3 ■ = 1 + = l,bj = 0); i.e., 
into itself, while, e.g., the operator o~ x , corresponding to (aj — 0,bj — 1), is mapped 
onto (aj = + 1 = 1, bj — 1), which coincides with o~ y . In the same way o~ y is mapped 
onto a x , while the identity (aj = 0, bj = 0) is mapped onto itself. 
In terms of field elements these transformations read 

z-rotation: « - « + ^tr^), 
P -> P, 
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Table 7. Curves and their corresponding transformations from a = 0, (3 = a n 
(left) and f3 = not (right). The x-, y and 2-rotations are indicated as x, y and z, 
respectively. 



Curve (ray) 



Transformation II Curve (ray) 



Transformation 



Q = ct 2 (k + re 2 ) 

ft = cr 2 K + K 2 

a = cr(re + re 2 ) 
- ere 2 



f} = C7K 
13 = 
13 = a 
(3 = era 



crf3 + f3 2 
cr 2 (3 + f3 2 



1 ® y 
j/® 1 

y®y 

z® z 
y®z 
z®y 
1 ® z 
z® 1 



a = re + re 
/3 = UK + o" 

/9 = cr 2 a 

13 = a 2 

f3 = a + a 2 

a = /3 + (3 2 



z(&y 



1 ® x 
1 ® 3/ 

y ® 1 



^ r ° tati0n: ^J+^trM,), (5 - 3) 

y-rotation- « - « + *i tr[(a + flty], 
y rotation. p ^ p + 6jt[ [((x + 0)9^ 



In particular, applying the above transformations to a ray (|4.1[) we get 
z— rotation 



(ri + (6 j )K + k 2 C 2 , 



[3^[3 = (k, 

or— rotation : „ f/ - , n \ i 2 (5-4) 



y— rotation 



/3->(C + C^+^> + « 2 (C + ?7) 2 , 



which are explicitly nonlinear operations. 

Note that the z- and ^-transformations produce regular curves starting from a 

ray 

z-rotation : a = i]C l P -> (vC 1 + 9j)P + P 2 , 

(5.5) 

x— rotation : /3 = f] Qa — > (i]^ 1 ^ + 6j)a + a 2 . 

Meanwhile, the y-rotation may lead to an exceptional curve (as it happens when we 
start with the horizontal or the vertical axes, £ = or r\ = 0). 

An important result to stress is that it is possible to obtain all the curves 
of the form (|4.6p and (|4.7p from the rays after some (nonlinear) operations (|5.4[) , 
corresponding to local transformations. The families of such transformations are the 
following: 

I. The rays and curves corresponding to factorable basis can be obtained from a 
single ray a — 0, = a 2 K (vertical axis) as shown in table 7 (left). 

II. The rays and curves corresponding to nonfactorable basis can be obtained 
from the ray a = an, (3 = ct 2 k (f3 = era) as shown in table 7 (right). 

This means that all the different tables can be generated from the standard 
one, given in table 1, by applying only local transformations that do not change 
the factorization properties of the MUBs. So, tables 2 to 6 are obtained from table 1 
from the transformations given in table 8. 
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Table 8. Transformation operators converting table 1 into each one of the tables 
indicated in the left column. Again, x-, y- and z-rotations are indicated as x, y 
and z, respectively. 



Table Transformation 



2 


X § 


9 1 


3 


1c 


5 z 


4 


y 5 


JZ 


5 


y 5 


5 x 


6 


1c 


%y 



The full set of striations for each bundle of curves (each table) is obtained by 
constructing "parallel curves" in the bundle in an obvious way: 

oi\{k) = ^ k + ^ik 2 , Px(k) = ijoK + rjiK 2 + A, (5.6) 
with A e GF(4). It is clear that no (a\(K) , 0\(k)) curve intersects the curve 
(a A '(/s), &'(*)) for A ^ A'. 

6. Extension to larger spaces 

The relation between Abelian curves in discrete phase space and different systems of 
MUBs can be extended to higher (power of prime) dimensions. For the most interesting 
n-qubit case, a generic Abelian curve (|4.2|) has the following parametric from 

n — 1 n — 1 

a{n) = ^ MmK 2 ™ , (3(k) = ^ rj m K 2m , (6.1) 

m— m— 

with fi m ,i] mi K £ GF(2 n ), and the commutativity condition takes now the invariant 
form 

5^tr(0m»fc)=O. (6.2) 

The simplest example of such curves are obviously the rays, parametrically defined 
as in Eq. (|4.1[) , where the conditions (|6.2[) arc trivially satisfied. Imposing the 
nonintersecting condition we can, in principle, get all the possible bundles of 
commutative curves. Nevertheless, in higher dimensions it is impossible to obtain 
all the curves from the rays by local transformations. This leads to the existence of 
different nontrivial bundles of nonintersecting curves, and consequently to MUBs with 
different types of factorization [22,36]. 

The problem of classification of bundles of mutually nonintersecting, nonsingular 
Abelian curves and its relation to the problem of MUBs in higher dimensions, and 
in particular the transformation relations between different MUB structures, will be 
considered elsewhere. 

7. Conclusions 

A new MUB construction has been worked out, with special emphasis in the two-qubit 
case. Its essential ingredient is a mapping between displacement operators, physical 
spin-1/2 operators and discrete phase-space curves. In phase space any nonsingular 
bundle of curves that fills every point and has only one common intersecting point 
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(here taken to be the origin) will map onto a MUB. The corresponding displacement 
operators can be obtained from these phase-space curves. 

For the two-qubit case, we have derived all the admissible curves and classified 
them into rays and curves (regular and exceptional, depending on degeneracy). In 
total, six different bundles can be constructed from the set of five rays and ten curves. 
We have also shown how the six tables representing sets of MUBs are related by local 
transformations, i.e., physical rotations around the x-, y- and z-axes. It is obvious 
that such rotations will not change the MUBs entanglement properties. 

A Wigner function can also be associated to each phase-space structure. Although 
we have not pursued this topic in the paper, it is straightforward to use any of the 
phase-space structures and follow the algorithm described in reference [33] (although 
in that paper the construction applies only to rays) to obtain such a function [39] . 

It is also formally straightforward to extend the method to any Hilbert space 
whose dimension is a power of a prime. However, only in the bipartite case one 
will find that all structures are related through local transformations. Already in the 
tripartite case different classes of entanglement exist [36] , and consequently some MUB 
structures are related through nonlocal (entangling) transformations. The extention 
of the present method provides a systematic way to find these transformations. 
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